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ABSTRACT 

Experiments designed to measure the redshifted 21 cm line from the Epoch of Reionization 
(EoR) are challenged by strong astrophysical foreground contamination, ionospheric distor- 
tions, complex instrumental response and other different types of noise (e.g. radio frequency 
interference). The astrophysical foregrounds are dominated by diffuse synchrotron emission 
from our Galaxy. Here we present a simulation of the Galactic emission used as a foreground 
module for the LOFAR- EoR key science project end-to-end simulations. The simulation 
produces total and polarized intensity over 10° x 10° maps of the Galactic synchrotron and 
free-free emission, including all observed characteristics of the emission: spatial fluctuations 
of amplitude and spectral index of the synchrotron emission, together with Faraday rotation 
effects. The importance of these simulations arise from the fact that the Galactic polarized 
emission could behave in a manner similar to the EoR signal along the frequency direction. 
As a consequence, an improper instrumental calibration will give rise to leakages of the po- 
larized to the total signal and mask the desired EoR signal. In this paper we address this for 
the first time through realistic simulations. 

Key words: cosmology: theory, diffuse radiation, observation; radio lines: general; 
instrumentation: interferometers; radio continuum: general 



1 INTRODUCTION 

The Epoch of Reionization (EoR) is expected to occur between 
redshift 6 and 12, as indicated from observed comic microwave 
background (CMB) polarization (Komatsu et al. 2009), high red- 
shift quasar spectra (Fan et al. 2006), and the thermal history of 
of the intergalactic medium (Theuns et al. 2002a,b; Bolton et al. 
2010). At redshifts, the 21 cm line from neutral hydrogen is shifted 
into meter wavelengths and therefore sets the frequency range of 
EoR experiments to the long- wavelength part of the radio spectrum 
(- 100 - 200 MHz). 

There are several planned and ongoing experiments designed 
to probe the EoR through redshifted 21 cm emission line from 
neutral hydrogen using radio arrays: GMRT\ LOFAR^, MWA^, 
21CMA^, PAPER^ and SKA^ 
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Giant Metrewave Telescope, http://gmrt.ncra.tifr.res.in 

Low Frequency Array, http://www.lofar.org 

Murchinson Widefield Array, http://www.mwatelescope.org/ 

21 Centimeter Array, http://21cma.bao.ac.cn/ 

Precision Array to Probe EoR, http://astro.berkeley.edu/"' dbacker/eor 

Square Kilometer Array, http://www.skatelescope.org/ 



The low-frequency radio sky at these wavelengths is domi- 
nated by diffuse synchrotron emission from the Galaxy and inte- 
grated emission from extragalactic sources (radio galaxies and clus- 
ters). Although, this foreground emission is 4 — 5 orders of magni- 
tude stronger than the expected EoR signal, the ratio between their 
intensity fluctuations on arcmin to degree scales measured by in- 
terferometers is 'only' 2 — 3 orders of magnitude (Shaver et al. 
1999). In addition to the foregrounds, the EoR experiments are also 
challenged by understanding of the instrumental response and iono- 
spheric disturbances to high precision (Labropoulos , in prepara- 
tion). 

Currently there are numerous efforts to simulate all the data 
components of the EoR experiments: cosmological 21 cm signal, 
foregrounds, ionosphere and instrumental response. The main aim 
of these end to end simulations is to develop a robust signal extrac- 
tion scheme for the extremely challenging EoR observations (e.g. 
Santos et al. 2005; Morales et al. 2006; Wang et al. 2006; Jelic et al. 
2008; Bowman et al. 2009; Harker et al. 2009a,b; Labropoulos et al. 
2009). 

The foregrounds in the context of the EoR measurements have 
been studied theoretically by various authors. Shaver et al. (1999) 
have given the first overview of the foreground components. Di 
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Matteo et al. (2002, 2004) have studied emission from unresolved 
extragalactic sources at low radio frequencies. Oh & Mack (2003) 
and Cooray (2004) have considered the effect of free-free emis- 
sion from extragalactic haloes. Santos et al. (2005) carried out a de- 
tailed study of the functional form of the foreground correlations. 
Jelic et al. (2008) have made the first detailed foreground model 
and have simulated the maps that include both the diffuse emis- 
sion from our Galaxy and extragalactic sources (radio galaxies and 
clusters). Gleser et al. (2008) have also studied both galactic and 
extragalactic foregrounds, de Oliveira-Costa et al. (2008) has used 
all publicly available total power radio surveys to obtain all- sky 
Galactic maps at the desired frequency range and Bowman et al. 
(2009) has studied foreground contamination in the context of the 
power spectrum estimation. 

Recently, a Galactic 3D emission model has been devel- 
oped by Sun et al. (2008); Waelkens et al. (2009); Sun & Reich 
(2009) (the Hammurabi'^ code), derived from a 3D distribution of 
the Galactic thermal electrons, cosmic-ray electrons and magnetic 
fields. The code is able to reproduce all-sky or zoom-in maps of the 
Galactic emission over a wide frequency range. 

In addition to simulations, a number of observational projects 
have given estimates of Galactic foregrounds in small selected 
areas. Ali et al. (2008) have used 153 MHz observations with 
GMRT to characterize the visibility correlation function of the fore- 
grounds. Rogers & Bowman (2008) have measured the spectral in- 
dex of the diffuse radio background between 100 and 200 MHz. 
Pen et al. (2009) have set an upper limit to the diffuse polar- 
ized Galactic emission; and Bernardi et al. (2009, 2010) obtained 
the most recent and comprehensive targeted observations with the 
Westerbork Synthesis Radio Telescope (WSRT). 

However, current observations are not able to fully constrain 
the foregrounds, especial the Galactic polarized synchrotron emis- 
sion, as required by EoR experiments. The importance of the polar- 
ized foreground stems from the fact that the LOFAR instrument, in 
common with all current interferometric EoR experiments, has an 
instrumentally polarized response. An improper polarization cali- 
bration will give rise to leakages of the complex polarized signal to 
the total signal. Since the Galactic polarized emission is quite struc- 
tured along the frequency direction, the leakage of polarized inten- 
sity will have similar structures along the frequency and will mimic 
the EoR signal. Therefore, for reliable detection of the EoR signal 
it is essential at this stage: (i) to simulate the polarized foregrounds, 
and (ii) to test the influence of the leakages on the extraction of the 
EoR signal. This paper focuses on both aspects. 

In our previous foreground model (Jelic et al. 2008), the total 
intensity Galactic emission maps were obtained from three Gaus- 
sian random fields. The first two were for the amplitudes of syn- 
chrotron and free-free emission and the third was for the spectral 
index of the synchrotron emission. The polarized maps were simu- 
lated in a similar way but with added multiple 2D Faraday screens 
along the line of sight. Despite the ability of that model to simulate 
observed characteristics of the Galactic emission (e.g. spatial and 
frequency variations of brightness temperature and its spectral in- 
dex), the model had some disadvantages: e.g. the Galactic emission 
was derived ad hoc and depolarization effects were not taken into 
account. 

Our new foreground model, presented in this paper, is an ex- 
tension of our previous foreground model (Jelic et al. 2008). It 
simulates both Galactic synchrotron and free-free emission in total 
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and polarized intensity, but in a more realistic way and with some 
additional benefits. The Galactic emission in our current model is 
derived from the physical quantities and 3D characteristics of the 
Galaxy (the cosmic ray and thermal electron density, and the mag- 
netic field; e.g. Beck et al. 1996; Berkhuijsen et al. 2006; Sun et al. 
2008, and references therein). In addition, the model has the flexi- 
bility to simulate any peculiar case of the Galactic emission includ- 
ing very complex 3D polarized structures produced by "Faraday 
screens" and depolarization due to Faraday thick layers. 

Our Galactic emission model has some similarities with the 
HAMMURABI model, but the difference between the two is the main 
purpose of the simulations. The HAMMURABI simulation is based 
on a very complex Galactic model with aim to reproduce the ob- 
served all-sky maps of the Galactic emission. Because of its com- 
plexity, the high resolution zoom-in maps require a lot of comput- 
ing power and time (Sun & Reich 2009). In contrast, our model 
is restricted to produce fast and relatively small maps of Galac- 
tic emission, which are then used as a foreground template for the 
LOFAR-EoR end to end simulation. Since the foreground subtrac- 
tion is usually done along the frequency direction, our model also 
includes 3D spatial variations of the spectral index of the Galactic 
synchrotron radiation. 

The paper is organized as follows. Section 2 gives a brief theo- 
retical overview of the Galactic emission and Faraday rotation. The 
observational constrains of the Galactic emission are presented in 
Sec. 3. The simulation algorithm is described in Sec. 4, while a 
few examples of simulated maps of the Galactic emission are pre- 
sented in Sec. 5. Section 5 also give a quantitative comparison of 
our simulations with the observations. Section 6 describes the EoR 
signal and instrumental response simulations of the LOFAR-EoR 
pipeline. We discuss the influence of the polarized foregrounds on 
the extraction of the EoR in Sec. 7. The paper concludes with sum- 
mary and conclusions (Sec. 8). 



2 THEORY 

In radio astronomy, at frequencies where the Rayleigh- Jeans law is 
applicable, the radiation intensity, / (energy emitted per unit time 
per solid angle and per unit area and unit frequency), at the fre- 
quency V is commonly expressed in terms of the brightness tem- 
perature (Tb): 



TbM 



r/M 



(1) 



where c is the speed of light and ks Boltzmann's constant. 

The emission coefficient, j (energy emitted per unit time per 
solid angle and per unit volume), at a certain frequency can also be 

2 

expressed in terms of the unit temperature, jh[T^) — 2k u'^ Ji^)' ^^ 
that: 



n{iy) 



/ ib(z^)ds, 
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where the integral is taken along the line of sight (LOS). 

In the following subsection we will give a brief theoretical 
overview of the Galactic synchrotron and free-free emission, as 
well as Faraday rotation, that will be used later in the simulation. 
The Galactic emission will be expressed in terms of jt and T^. 

2.1 Synchrotron emission 

Synchrotron emission originates from the interaction between rela- 
tivistically moving charges and magnetic fields. In our own galaxy. 
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synchrotron emission arises from cosmic ray (CR) electrons pro- 
duced mostly by supernova explosions and the Galactic magnetic 
field. A fairly complete exposition of the synchrotron emission the- 
ory is presented in e.g. Pacholczyk (1970) and Rybicki & Lightman 
(1986). Here we only give a simple description of the emission. 

The Galactic synchrotron emission is partially linearly polar- 
ized. Its properties depend on the spatial and energy distribution 
of the CR electrons, and the strength and orientation of the per- 
pendicular (with respect to the LOS) component of the Galactic 
magnetic field, B± . The emission coefficients of the Galactic total 
and polarized synchrotron radiation, jl^^^ and j^^*^^, are given 
respectively in cgs units, at the frequency zy, by: 
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The charge of the electron is given bye = 4.810~^° Fr, the mass 
by rUe = 9.1 • 10~^^ g and ucr is the CR electron density. Note 
that for the CR electrons we assume that their energy spectrum is a 
power law with a spectral index p: N{j)d'j = ^CRo7~^» where 7 
is the Lorenz factor and A^(7) the number density of electrons with 
energy between 7 and 7 + d7 and ncRo normalization constant. 
Furthermore, we assume that their velocity and pitch angle distri- 
bution is isotropic. Both simplifications are consistent with obser- 
vations and are widely used by many authors (e.g. Sun et al. 2008; 
Waelkens et al. 2009, as most recent examples). Note that the in- 
trinsic degree of polarization of synchrotron radiation depends on 
the energy spectral index p and is given by 
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p+1 



(6) 



p + 7/3' 

The Stokes Q and U parameters of the polarized Galactic syn- 
chrotron emission are given by: 

(7) 
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where $ is polarization angle defined with respect to the orientation 
of the magnetic field. 

By integrating j^^ ^^^, j^ & j^ along some LOS (see Eq. 2) we 
get the total and polarized Galactic synchrotron emission in terms 
of the brightness temperature (T/^^"", T^ & T^). Note that ob- 
served polarized emission and polarization angle ^obs are given 
by: 
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2.2 Free-free emission 
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Radiation due to (de) acceleration of a charged particle in the elec- 
trical field of another is called bremsstrahlung or free-free radia- 
tion. The Galactic free-free emission originates from electron-ion 
encounter in the warm ionized gas. As for the synchrotron emis- 
sion a detail theory of the free-free radiation can be found in e.g. 



Rybicki & Lightman (1986) and Wilson et al. (2009), here we give 
only the necessary formulae. 

The optical depth, r^^, of the warm ionized gas at a given low 
radio frequency v is: 
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where Te is temperature of the ionized gas, g/f is the Gaunt factor 
of the free-free transition given by: 
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and EM is emission measure defined as: 
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The integral is taken over the LOS, where Ue in cm~^ is the elec- 
tron density of the warm ionized gas. 

The Galactic free-free emission in terms of brightness temper- 
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ature, jf^-' , is given by: 
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Note that for optically thin ionized gas, j^^ is simply given by 



2.3 Faraday rotation 

When the polarization angle of an electromagnetic wave is rotated 
while passing through a magnetized plasma, the effect is called 
Faraday rotation (for details see Rybicki & Lightman 1986; Wilson 
et al. 2009). The rotation depends on the frequency of the wave, u, 
electron density, rie, and magnetic field component parallel to the 
LOS, B,,: 
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where the polarization angle of the wave before rotation is denoted 
with $0- Eq. 15 is also written as $ = $o + RMX^ with A in units 
of m and RM (rotation measure) defined as: 

B\\ ds 
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The RM is positive when B\\ points towards observer and negative 
when B\\ points in away. 



3 OBSERVATIONAL CONSTRAINTS 

There are several all-sky maps of the total Galactic diffuse radio 
emission at different frequencies and angular resolutions (Haslam 
et al. 1982; Reich & Reich 1986, 1988; Page et al. 2007). The 
150 MHz map by Landecker & Wielebinski (1970) is the only all- 
sky map in the frequency range (100 — 200 MHz) relevant for the 
EoR experiments, but has only 5° resolution. 

At high Galactic latitudes the minimum brightness tempera- 
ture of the Galactic diffuse emission is about 20 K at 325 MHz 
with variations of the order of 2 per cent on scales from 5 to 30 
arcmin across the sky (de Bruyn et al. 1998). At the same Galactic 
latitudes, the temperature spectral index of the Galactic emission 
is about -2.55 at between 100 and 200 MHz (Rogers & Bow- 
man 2008) and steepens towards higher frequencies (e.g. Platania 
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et al. 1998; Bennett et al. 2003; Bernardi et al. 2004). Furthermore, 
the spectral index gradually changes with position on the sky. This 
change appears to be caused by a variation in the spectral index 
along the line of sight. An appropriate standard deviation in the 
power law index, in the frequency range 100-200 MHz appears to 
be of the order of - 0.1 (Shaver et al. 1999). 

Using the obtained values at 325 MHz and assuming the fre- 
quency power law dependence, the Galactic diffuse emission is ex- 
pected to be 140 K at 150 MHz, with -- 3 K fluctuations. 

Studies of the Galactic polarized diffuse emission are done 
mostly at high radio {^^ IGHz) frequencies (for a recent re- 
view see, Reich 2006). At lower frequencies {r^ 350 MHz), there 
are several fields done with the Westerbork telescope (WSRT) 
(Wieringa et al. 1993; Haverkorn et al. 2003; Schnitzler 2008). 
These studies revealed a large number of unusually shaped polar- 
ized small-scale structures of the Galactic emission, which have no 
counterpart in the total intensity. These structures are usually at- 
tributed to the Faraday rotation effects along the line of sight. 

At high Galactic latitudes, the Galactic polarized emission at 
350 MHz is around 5 K or more, on 5-10 arcmin scales (de Bruyn 
et al. 2006). At 150 MHz this polarized emission would scale to 
few tens of Kelvin if it were Faraday thin. However, depolarization, 
that is prominent at low radio frequencies, can significantly lower 
the level of polarized emission. 

Recently, a comprehensive program was initiated by the 
LOFAR-EoR collaboration to directly measure the properties of 
the Galactic radio emission in the frequency range relevant for the 
EoR experiments. The observations were carried out using the Low 
Frequency Front Ends (LFFE) on the WSRT radio telescope. Three 
different fields were observed. The first field was a highly polarized 
region known as the "Fan region" in the 2nd Galactic quadrant at a 
low Galactic latitude of '-^ 10° Bernardi et al. (2009). The second 
field was a very cold region in the Galactic halo (/ ^ 170°) around 
the bright radio quasar 3C196, and third was a region around the 
North Celestial Pole (NCP, / - 125° Bernardi et al. 2010). The last 
two fields represent possible targets for the LOFAR-EoR observa- 
tions. Below we present the main results of these papers. 

In the "Fan region", fluctuations of the Galactic diffuse emis- 
sion were detected at 150 MHz for the first time. The fluctuations 
were detected both in total and polarized intensity, with an rms of 
14 K (13 arcmin resolution) and 7.2 K (4 arcmin resolution) re- 
spectively (Bernardi et al. 2009). Their spatial structure appeared 
to have a power law behavior with a slope of —2.2 ± 0.3 in total in- 
tensity and — 1 .65 =b 0. 15 in polarized intensity. Note that, due to its 
strong polarized emission, the "Fan region" is not a representative 
part of the high Galactic latitude sky. 

Fluctuations of the total intensity Galactic diffuse emission in 
the "3C196" and "NGP" fields were also observed on scales larger 
than 30 arcmin, with an rms of 3.3 K and 5.5 K respectively. 

Patchy polarized emission was found in the "3C196" field, 
with an rms value of 0.68 K on scales larger than 30 arcmin 
(Bernardi et al. 2010). Thus, the Galactic polarized emission fluc- 
tuations seem to be smaller than expected by extrapolating from 
higher frequency observations. Recent observations at mid-galactic 
latitude with the Giant Metrewave Radio Telescope (GMRT) con- 
firm this conclusion, by setting an upper limit to the diffuse polar- 
ized Galactic emission in their field to be < 3 K at 150 MHz and 
on scales between 36 and 10 arcmin (Pen et al. 2009). 



4 SIMULATION 

In this section, the various components of the simulation that 
lead towards the brightness temperature maps of the Galactic syn- 
chrotron and free-free emission in a total and polarized intensity 
are explained. Because the simulated maps will be used as a fore- 
ground template for the LOFAR-EoR end-to-end simulations, the 
foreground simulations assume the angular and frequency range of 
the LOFAR-EoR experiment, i.e. 10° x 10° maps from 115 MHz 
to 180 MHz. In addition, all parameters of the simulations can be 
tuned to any desired value or have any desired characteristic, allow- 
ing to explore the parameter space of our Galactic model. 

Our algorithm is based on a 3D grid in a Cartesian coordi- 
nate system, where x^-plane represents the angular plane of the 
sky ("flat sky" approximation valid for a small field of view) and z 
axis is a line of sight direction in parsecs. 

The first step in our simulation is to calculate, at a certain fre- 
quency, the 3D emission coefficient of the Galactic synchrotron and 
free-free emission expressed in terms of unit temperature (see Eq. 3 
& 14). The emission coefficients are obtained from the cosmic-ray, 
ncR, and thermal electron, rie, densities, and the Galactic mag- 
netic field {B). Given the 3D emission coefficients, we integrate 
along the LOS to obtain the brightness temperature maps of the 
Galactic synchrotron and free-free emission at a certain frequency. 
The calculation also includes Faraday rotation effects. Note that all 
parameters of the simulation are set in such a way that simulated 
maps are quantitatively (e.g. presence of the structures at different 
scales, spatial and frequency variations of the brightness temper- 
ature and its spectral index, etc.) in agreement with the observa- 
tions overviewed in Sec. 3. A detail description of all input param- 
eters (ncR, rie and B) and the algorithm follow. A few examples 
of simulated Galactic emission data cubes are presented in Sec. 5, 
together with a quantitative comparison with the observations. 



4.1 Cosmic ray electron density 

The cosmic ray (CR) electrons relevant for the Galactic synchrotron 
emission have energies between 400 MeV and 25 GeV, assuming 
a Galactic magnetic field of a few fiG (Webber et al. 1980). In 
this energy range, the CR electron distribution can be described as 
a power law. The power law is normalized according to the mea- 
surements obtained in the solar neighborhood. However, the locally 
measured values might not be a good representative for the CR den- 
sity elsewhere in the Galaxy (e.g. Strong et al. 2004). As a conse- 
quence, the CR electron distribution is weakly constrained. 

In our simulation, uniform CR electron density distribution is 
assumed in the x2/-plane. In the z direction we follow Sun et al. 
(2008) and assume an exponential distribution: 



^CR = ^CRo exp 



VlkpcJ- 



(17) 



Note that ncRo depends on the assumed energy spectral index p 
of the CR electrons, so it is normalized according to Eq. 2 for the 
synchrotron radiation. Assuming Tb(150 MHz) :^ 145 K, B± — 
5 /xG andp = 2, we get ncRo — 1-4 • 10~^ cm~^. 

In the desired frequency range of our simulation, the assumed 
energy spectral index p = 2 is consistent with the values of the typ- 
ically observed brightness temperature spectral index of the Galac- 
tic synchrotron emission {(3 — —2.5, see Sec. 3 )^. 



The brightness temperature spectral index (3 of the Galactic synchrotron 
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In addition, the spatial variations of the spectral index p are in- 
troduced to mimic the observed spatial fluctuations of /3. We follow 
our previous model (Jelic et al. 2008) and simulate the variation of 
p (or 13) as a Gaussian random field. For the power spectrum of the 
Gaussian random field it is assumed a power law with index —2.7 
(JeUc et al. 2008). 



4.2 Galactic magnetic field 

The Galactic magnetic field has two components: a regular com- 
ponent Br and a random component h, so that the total Galactic 
magnetic field is given as 5 = Br + 6 (for review see Beck et al. 
1996; Han & Wielebinski 2002). The regular component is usually 
simulated as a combination of a disk and a halo field, whereas the 
random field component is simulated as a Gaussian random field, 
Gaussian random field, (for details see Sun et al. 2008; Sun & Re- 
ich 2009). Note that for our calculations, we split 5 in a component 
parallel {B\\) and perpendicular {B±_) to the LOS, so that Faraday 
rotation is defined by B\\ and synchrotron emission by B±. 

Considering the aim of our effort to simulate the Galactic 
emission for a small patch of the sky, we treat the regular field com- 
ponent in a simplified way. The regular field component is assumed 
to be uniform in the x^y-plane and to have an exponential decrease 
in the z direction. The typical value of the regular field component 
is a few fiG (for review see Beck et al. 1996; Han & Wielebinski 
2002). 

For the random field component we follow Sun et al. (2008); 
Sun & Reich (2009) and simulate it as a Gaussian random field. The 
power spectrum of the field follows a power law, with spectral index 
—8/3. This spectral index is commonly used for a Kolmogorov- 
like turbulence spectrum. 

In our simulation, the realization of the random field compo- 
nent is done in the following way. First we generate three differ- 
ent Gaussian random fields for the bx, by and bz component. From 
those three fields, we then calculate the amplitude of b and normal- 
ize it to the desired value. A typical value for the mean random field 
strength is 6 = 3 /iG (Sun et al. 2008). 



4.3 Thermal electron density 

At high Galactic latitudes, the warm ionized medium consists 
mostly of diffuse ionized gas (DIG) with total emission measure of 
- 5 pc cm-^ and Te = 8000 K (Reynolds 1990). The properties 
of the DIG can be traced by its free-free emission and dispersion 
measure^ (DM) of pulsars (e.g. Gaensler et al. 2008). 

Recent simulations of the Galactic emission (Sun et al. 2008; 
Waelkens et al. 2009; Sun & Reich 2009) used the Cordes & Lazio 
(2002) model for the thermal electron distribution. That model sim- 
ulates the Galaxy as several large-scale (e.g. thin and thick disk, and 
spiral arms) and small-scale (e.g. supernovae bubbles) structures. 
In our simulation, we follow our previous model of the Galactic 
free-free emission (Jelic et al. 2008) and simulate the thermal elec- 
tron density distribution as a Gaussian random field with the power 
law type of the spectrum. The spectral index of the power law is 



emission and the energy spectral index p of the CR electrons are related as 
l3 = -{p + S)/2. 

^ The dispersion measure is defined as the integral of the thermal electron 
density along the LOS. Knowing the distance to the pulsars (e.g. determined 
by parallax), an electron density model can be obtained by a fit to the ob- 
served DMs. 
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Figure 1. Flow chart of the algorithm: the Galactic emission (synchrotron 
and free- free) is derived from the physical quantities and 3D characteris- 
tics of the Galaxy, i.e. cosmic ray, ncR, and thermal electron, Ue, density; 
and magnetic field, B. In addition, the algorithm includes Faraday rotation 
effects. 



—3. The amplitude of the Gaussian random field (thermal electron 
density) is normalized in a way to match the typical observed EM 
of the quasars at high Galactic latitudes (EM values are taken from 
Berkhuijsen et al. 2006). 

It is important to note that our model is flexible to include ad- 
ditional features of the thermal electron distribution, e.g. dense bub- 
bles or clumpy distribution. Some of these features are presented in 
Sec. 5. 



4.4 The Algorithm 

Here we summarize the steps we follow to obtain maps of the 
Galactic emission at a desired frequency. The flow chart of the al- 
gorithm is presented in Fig. 1. 

(i) The CR electron density, ncR, and the regular component, 
Br, of the Galactic magnetic field are defined on 3D grid. The dis- 
tributions of ncR and Br are uniform in the xy-p\miQ and have an 
exponential decrease in the z direction. 

(ii) The spatial distribution of the CR electron energy spectral 
index, p, the random component, b, of the Galactic magnetic field 
and the thermal electron density, rie, are simulated as Gaussian ran- 
dom fields. The Gaussian random fields are normalized to result in 
a desired rms value of the brightness temperature maps. Note that 
additional features in the electron distribution are added if desired. 

(iii) The parallel, B\\ , and perpendicular, B±, component of the 
total Galactic magnetic field, B, are calculated from Br and b. 

(iv) Using Eq. 3, the emission coefficients of the Galactic total, 
j^ *^^, and polarized, j^ *^^, synchrotron radiation are calculated. 

(v) The optical depth, r^^, and emission coefficient, j^^ , of the 
thermal plasma are obtained from Eq. 11 & 14. Note that these 
effect is really significant only on the lowest radio frequencies. 

(vi) Absorption of the synchrotron emission by the optical thick- 
ness of the ionized plasma is taken into account as exp {—r^^) 
factor. 

(vii) Using Eq. 15, the Faraday rotation effect is calculated and 
the polarization angle, $, is obtained. Note that the intrinsic polar- 
ization angle, $o, is defined as the inclination of B±. 

(viii) The Stokes Q and U emission coefficients of the polarized 
emission, j^ & j^, are calculated using Eq. 1 &^. 

(ix) By integrating jl^^^,j^ & jb along some LOS (see Eq. 2), 
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the total and polarized Galactic synchrotron emission in terms of 
the brightness temperature, t/^^"", T^ Sc T^ , are obtained. 

(x) Finally, the maps of the total polarized emission (Tf^) and 
observed polarization angle ^obs is calculated using Eq. 9 & 10. 

In the following section we will show some examples of the Galac- 
tic emission maps obtained by this algorithm. 



5 EXAMPLES OF SIMULATED GALACTIC EMISSION 

Here we demonstrate the ability of our algorithm to realistically 
simulate Galactic synchrotron and free-free emission both in total 
and polarized intensity. We also quantitatively compare results of 
our simulations with the observations. 
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5.1 Simulated data cubes 

The simulated data cubes are presented for four cases of the Galac- 
tic emission. First three models show ability of our algorithm to 
simulate different cases of Faraday rotation and depolarization, 
while the last model is tailored to be in agreement with Bemardi 
et al. (2009) observations of the Fan region. 

• Model A: CR electrons are distributed in a region of 1 kpc in 
depth along the LOS (ncRo = 1.4 x 10"^ cm~^). In front of this 
region, there is a thermal electron cloud of 300 pc in depth along 
the LOS, with an average emission measure of EM = 8 cm~^pc, 
and Te = 8000 K. The thermal electron cloud is acting as a "Fara- 
day screen" that rotates the polarization angle of synchrotron emis- 
sion. It is assumed that Galactic magnetic field is uniform through- 
out the region (Br,\\ = 3 /xG and Br,± = 2.5 /xG). 

• Model B: Both CR and thermal electrons are mixed in the re- 
gion of 1 kpc in depth along the LOS (ncRo, EM, Te and Br 
are the same as in model A). The polarized synchrotron radiation is 
differentially Faraday rotated. Since the Galactic magnetic field has 
the same characteristics as in model A, depolarization is produced 
only by thermal plasma. 

• Model C: The same as model B, but in the middle of the sim- 
ulated region there is a dense thermal electron bubble 300 pc in 
depth along the LOS, with a strong magnetic field (B\\ = 10 /xG). 
The Faraday rotation along the bubble will be much larger than 
in other parts of the region. Note that the size of the bubble is 
quite large in order to make its appearance in the final maps more 
clear. Moreover, it is assumed that Galactic magnetic field in not 
uniform throughout the region: Br,\\ exponentially decreases along 
the LOS. We also include a random magnetic field component as a 
Gaussian random field (b\\ = 3 /iG, see Sec. 4). Depolarization is 
then produced both by thermal plasma and variations of the mag- 
netic field. 

• Model D: This model is based on a proposed cartoon of the 
magnetized interstellar medium in the Fan region by Bernardi et al. 
(2009, see fig. 12). The intrinsically polarized background emis- 
sion passes through multiple Faraday screens which make the spa- 
tial distribution of the polarized emission more clumpy. A bubble 
of ionized thermal plasma sits between the Faraday screens and the 
observer. This bubble further rotates the polarized emission to the 
higher values of RM. All parameters of this model are set to match 
the observed properties of this region: fluctuations of the diffuse 
emission at 150 MHz and their power-law type of the power spec- 
trum (see Bernardi et al. 2009, and Sec. 3). 



Figure 2. The intensity distributions (see Fig. 3 & 4) are obtained for four 
different models of the Galactic emission . The first (model A) assumes that 
synchrotron and free-free emitters are spatially separated, so that thermal 
plasma acts as a "Faraday screen". The second, third and fourth (model B, 
C, & D) have regions where both types permeate in a different way. The 
synchrotron emission is differentially Faraday rotated and depolarization 
occurs. 



Cartoons of these four models are given in Fig. 2. The Galactic 
emission maps are obtained in the frequency range from 115 to 
180 MHz, with 0.5 MHz step. 

In the first three models the synchrotron emission originates 
from the same CR electron distribution with uniform B± com- 
ponent of the Galactic magnetic field. Therefore, the brightness 
temperature maps of the total and intrinsic ^° polarized synchrotron 
emission are equivalent in all three models. The same is valid for 
the free-free emission, i.e., rie is normalized to the same value of 
EM. The resulting 10° x 10° maps at 150 MHz are shown in 
Fig. 3.a, b & d, while the mean and rms of the maps are given in 
Table 1. A map of the brightness temperature spectral index /3 of 
simulated total intensity synchrotron emission is shown in Fig. 3.c. 

Simulated polarized emission maps of the first three Galactic 
synchrotron emission models are shown in Fig. 4. a, b & c. Their 
mean and rms values together with the degree of polarization and 
depolarization are listed in Table 1. A random line through syn- 
chrotron total and polarized intensity frequency data cubes of first 
three models are given in Fig. 6. 

Simulated total and polarized emission images of the Fan re- 
gion (model D) are given at 150 MHz in Fig. 5 a & b. For easier 
comparison with the observation an RM image, obtained by apply- 
ing rotation measure synthesis technique (Brentjens & de Bruyn 
2005) on the simulated data, is shown in Fig. 5 c. 



5.2 Discussion and tests of the algorithm 

First test of the algorithm is to estimate the degree of intrinsic po- 
larization, H, of the simulated map (Fig. 3.b). By dividing the 
intrinsic polarized emission map and the total intensity map, we 



^^ Here, the intrinsic polarized emission, iPIsyn, means emission defined 
by Eq. 3. Note that the polarization angle of this emission is assumed to be 
uniform across the whole region. Thus any effect caused by thermal elec- 
trons or Bu will be immediately apparent. 
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obtain 11 = 0.69 (see Table 1). This value is in a good agreement 
with the expected theoretical value 11^=2 = 9/13 (see Eq. 6). 

Further we explore morphological characteristics of the sim- 
ulated polarized emission by comparing the image of intrinsically 
polarized emission (Fig. 3.b) with the images of polarized emission 
of the first three models (see Fig. 4). Note that degree of polarized 
and depolarized emission is given in Table 1. 

Model A assumes that there is no region in which the plasma 
(thermal electron cloud) is mixed with CR electrons. Therefore, po- 
larization angle of synchrotron emission along the LOS are Faraday 
rotated by an equal amount (defined by Eq. 15). Since there is no 
differential Faraday rotation, depolarization does not occur and the 
polarized intensity of the synchrotron emission is unchanged (see 
Fig. 4. a and Table 1). Morphology of the polarization angles over 
the map is determined only by the spatial RM variations of the 
plasma. 

In order to test the above results, we correlate a polarized 
emission image with a total intensity, and polarization angles with 
the image of free-free emission. The cross correlation coefficient^^ 
is C = 1 in both cases. Since the polarized emission s completely 
correlated with the total intensity, and the free-free emission, which 
traces thermal plasma, with polarization angles, we conclude that 
our algorithm is preforming as expected. The spatial structures of 
polarized emission are unchanged, while the structures of polariza- 
tion angle follow thermal plasma. 

In B and C models there are regions where both CR and ther- 
mal electrons are mixed. In those regions the polarization angle 
of the synchrotron radiation is then differentially Faraday rotated 
along the LOS. As a result, the polarized synchrotron radiation is 
quenched and the level of the polarized emission is weaker than in- 
trinsic polarization emission (see Table 1). In addition to depolar- 
ization produced by the mixed regions, in model C depolarization 
is also produced by variations of the Galactic magnetic field along 
the LOS. 

As for model A, we calculate correlations between the to- 
tal and polarized images. In both models (B & C) we find much 
weaker correlation (Cb,c = 0.3,0.1) than in model A. The 
weaker correlation results from differential Faraday rotation along 
the LOS, which changes the morphology of polarized emission 
both in intensity and angle. To illustrate this we correlate the po- 
larized synchrotron emission with the free-free emission. We find 
anti-correlation in both models (Cb,c — —0.7,-0.5). The inten- 
sity of the polarized emission is more quenched in the regions that 
show stronger free-free emission (regions of denser plasma). More- 
over, model C shows weaker (anti)correlations than model B, due 
to Faraday structures produced by magnetic field. 

The last aspect of the algorithm, which needs to be verified, 
is the frequency behavior of the simulated Galactic emission. It is 
expected that the synchrotron emission shows roughly a power law 
along the frequency (see Eq. 3). By examining Fig. 6, one can see 
that simulated synchrotron emission indeed shows a power law be- 
havior in total intensity. However, in polarized intensity only model 
A shows a power law behavior. All the other models show charac- 



^^ The cross correlation coefficient between two images {ai^j and hi^j) 
with the same total number of pixels n is defined as: 
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Table 1. The mean and rms value of the maps shown in the Fig. 3 & 4. 
All the values are given in kelvin. For completeness, degree of polarized 
{PI /T) and depolarized (dep.) emission is calculated. 
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teristic fluctuations caused by differential Faraday rotation (Bum 
1966). Note that fluctuations show richer structures, as in model C, 
if the differential Faraday rotation is more prominent. 

The calculated brightness temperature spectral index, Psyn, 
variations of the simulated Galactic synchrotron emission are given 
in the Fig. 3.c. The variations of ^syn are estimated by fitting the 
data with a linear functions along the frequency in logarithmic 
scale. The mean value of Psyn is -2.50, which is in a good agree- 
ment with the expected theoretical value /3p=2 = — (p + 3)/2 = 
—2.55. A slight difference between the two is caused by the 3D 
spatial variations of the spectral index p. 

Based on discussion and tests of the algorithm presented in 
this subsection, we can conclude that our algorithm is simulating 
Galactic emission as expected. In the following subsection, we will 
compare the simulated data cubes with the observations. 



5.3 Comparison with observations 

Input parameters of the simulation (properties of the CR electrons, 
thermal electrons, and Galactic magnetic field, see subsections 4.1, 
4.2 & 4.3) are chosen in such way that the global properties (bound- 
ary conditions) of the simulated maps, e.g. morphology and fre- 
quency behavior of the Galactic emission, are in agreement with 
the observations overviewed in Sec. 3: 

• the mean value of the simulated synchrotron emission in to- 
tal intensity, Isyn = 142 K (see Table 1) is comparable to ob- 
served emission at 150 MHz (in the cold regions of the Galaxy at 
the high galactic latitudes, Landecker & Wielebinski 1970) or to 
the values obtained by extrapolation from the higher frequencies, 
i.e.,- 150 K; 

• the mean and variations of the brightness temperature spec- 
tral index of the simulated maps, Psyn = —2.50 (see Subsec- 
tion 5.2) and a — 0.1, are in agreement with the observed values. 






-2.55 and CT? 



: 0.1 (e.g. Rogers & Bowman 2008, and 



where a (b) is the mean and a a (crb) the standard deviation of the image a 

(by 



references therein); 

• the spatial structures of the simulated Galactic emission in 
total intensity are morphologically similar to the observed emis- 
sion around 150 MHz (Bernardi et al. 2009, 2010) and - 300 - 
400 MHz (e.g. Haslam et al. 1982; de Bruyn et al. 1998, and refer- 
ences in Sec. 3), i.e., power law type of the power spectrum with a 
negative power law index (more power on the large scales); 

• the simulated polarized emission (e.g. model B, C & D) shows 
small-scale structures, which have no counter parts in the total in- 
tensity, e.g. like a number of observations with the WSRT (see 
Sec. 3); 

• the simulated distribution of the thermal electrons reproduce 
the typical observed EM ^ 10 cm~^pc at the high Galactic lati- 
tudes (Berkhuijsen et al. 2006). 

• the simulated Galactic emission (see Fig. 6) shows charac- 
teristic power law behavior along the frequency in total intensity. 
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(a) 



(b) 



(c) 



(d) 



Figure 3. Simulated maps of the total (Fig. 3. a) and intrinsic polarized (Fig. 3.b) intensity of the Galactic synchrotron emission. The polarization angle is 
plotted over the polarized map as a white lines. The map of brightness temperature spectral index (5 of simulated total intensity synchrotron emission is shown 
in Fig. 3.C. The total intensity map of the free-free emission is shown in Fig. 3.d. The angular size of the maps are 10° x 10°, with ~ 1 arcmin resolution. 
The color bar represents the brightness temperature T^ of emission in kelvin at 150 MHz. The mean and rms value of the maps are given in the Table 1. 





(a) 



(b) 



(c) 



Figure 4. Simulated maps of the polarized intensity and polarization angle (white lines) of the four different Galactic synchrotron emission models (A, B and 
C from left to right). The angular size of the maps are 10° x 10°, with ~ 1 arcmin resolution. The mean and rms values of the maps at 150 MHz are 
given in the Table 1 . 



while in polarized intensity can reproduce the characteristic polar- 
ization functions (Burn 1966). 

We also quantitatively compare the simulated emission in the 
Fan region (our model D) with the observations of this region at 
150 MHz (Bernardi et al. 2009): 

• simulated fluctuations have an rms of 14 K in total intensity 
and rms of 7 K in polarized intensity (these values are the same as 
the observed ones); 

• the power law indices of the power spectrum obtained from 
the simulation are in agreement with the observations: —2.2 in total 
intensity, and —1.6 in polarized intensity; 

• simulated maps (see Fig. 5) show similar morphological struc- 
tures (their statistical distribution of power) as the observed ones 
(fig 5., 8. & 10. in Bernardi et al. 2009); 

• a simulated rotation measure cube^^ shows similar morpho- 



^'^ The rotation measure cube is obtained using the rotation measure syn- 
thesis method (Brentjens & de Bruyn 2005). 



logical structures as the observed ones, i.e., at —2 rad m ^ there 
is a hole in the emission, a bubble with a diameter of ^ 2°. 



Here we emphasize that we have also successfully tested our 
algorithm at the higher frequencies (^ 350 MHz) by simulating 
Galactic emission towards the cluster Abel 2255. The simulated 
maps have been compared with the observations obtained by R.F. 
Pizzo et al. {in prep., private communication). Details of these sim- 
ulations will be presented in a separate paper (Jelic et al., in prep.). 

To summarize, simulated maps show all observed character- 
istics of the Galactic emission, e.g. presence of the structures at 
different scales, spatial and frequency variations of the brightness 
temperature and its spectral index, complex Faraday structures, and 
depolarization. Based on these results, we conclude that our model 
is able to simulate realistic maps of the Galactic emission both in 
total and polarized intensity, and can be used as an realistic fore- 
ground model in simulations of the EoR experiments. 
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Figure 5. Simulated maps of the Fan region (model D) in total (Fig. 5 a) and polarized intensity (Fig. 5 b). An RM image obtained by applying rotation 
measure synthesis technique on the simulated data is shown in Fig. 5 c. 
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Figure 6. A random line through synchrotron total intensity (Tlsyn) and polarized intensity (PIsyn, Qsyn) frequency data cubes. The solid lines are for the 
Galactic model A, dashed lines for model B and dotted lines for model D. Tlsyn is the same in all four models, while Usyn is not presented since it is 
similar to Qsyn- Note the polarized structures along the frequency direction. An improper polarization calibration of the instrument could cause a leakage of 
these structures to the total intensity and severely contaminate the EoR signal. 



6 LOFAR-EOR SIMULATION PIPELINE 



et al. 2008) and the instrumental response (described on Labropou- 
los et al. 2009). 



The LOFAR-EoR project relies on a detailed understanding of as- 
trophysical and non-astrophysical contaminations that can contam- 
inate the EoR signal: the Galactic and extragalactic foregrounds, 
ionosphere, instrumental effects and systematics. In order to study 
these components and their influence on the detection of the 
EoR signal, a LOFAR-EoR simulation pipeline is being devel- 
oped by the LOFAR-EoR team. The pipeline consists of the three 
main modules: the EoR signal (based on simulations described in 
Thomas et al. 2009), the foregrounds (based on this paper and Jelic 



In this paper we use the LOFAR-EoR pipeline to illustrate the 
need for excellent calibration of the instrument in order to reliably 
detect the EoR signal (see Sec. 7). The following two subsections 
present a brief overview of the EoR signal and instrumental re- 
sponse modules. 
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6.1 EoR signal 

The predicted differential brightness temperature of the cosmologi- 
cal 21 cm signal with the CMB as the background is given by (Field 
1958, 1959; Ciardi & Madau 2003): 



5Tu 



26mKxHi(l + ^) 1 



0.02 



10 



0.3 



1/2 



(19) 



Here Tg is the spin temperature, xhi is the neutral hydro- 
gen fraction, 5 is the matter density contrast and h — 
ifo/(100 kms~^Mpc~^). Throughout we assume ACDM- 
cosmology with WMAP3 parameters (Spergel et al. 2007): h — 
0.73, Qb = 0.0418, ^m = 0.238 and Qa = 0.762. In addition we 
assume that Ts ^ Tcmb, which is assumed in most of the current 
simulations. 

The cosmological 21 cm maps {5Th) are simulated using the 
BEARS algorithm (Thomas et al. 2009). BEARS is a fast algorithm 
to simulate the underlying cosmological 21 cm signal from the 
EoR. It is implemented using an N-body/SPH simulation in con- 
junction with a 1-D radiative transfer code under the assumption 
of spherical symmetry of the ionized bubbles. The basic steps of 
the algorithm are as follows: first, a catalogue of ID ionization 
profiles of all atomic hydrogen and helium species and the tem- 
perature profile that surround the source is calculated for different 
types of ionizing sources with varying masses, luminosities at dif- 
ferent redshifts. Subsequently, photon rates emanating from dark 
matter haloes, identified in the N-body simulation, are calculated 
semi-analytically. Finally, given the spectrum, luminosity and the 
density around the source, a spherical ionization bubble is embed- 
ded around the source, whose radial profile is selected from the 
catalogue as generated above. For more details we refer to Thomas 
et al. (2009). 

For the purpose of this paper we use the 5Tb data cube (2D 
slices along the frequency/redshift direction) of the cosmological 
21 cm signal for the 'Stars' patchy reionization model (see Thomas 
et al. 2009). The data cube consists of 850 slices in the frequency 
range from 115 MHz to 200 MHz with 0.1 MHz step (correspond- 
ing to redshift between 6 and 11.5). Slices have a size of 100 h~^ 
comoving Mpc and are defined on a 512^ grid. An example of a 
random line of sight through simulated 21 cm data cube, with an- 
gular and frequency resolution matching that of LOFAR, is shown 
in Fig. 7. 



6.2 Instrumental response 

In order to produce the dirty maps of the diffuse emission, we cal- 
culate the 2D Fourier transform of the data for each correlation on 
a fine grid of 1.2 arcmin. We assume that there are 24 stations in 
the LOFAR array that are used for the observations. We then use a 
bilinear interpolation to estimate the values of the visibilities at the 
uvw points that correspond to the points sampled by the interfer- 
ometer pairs of the core. This is done for 4hrs of synthesis, 10 sec 
integration and for the whole frequency range between 115 MHz 
and 180 MHz, using a step of 0.5 MHz. The above procedure is im- 
plemented as a parallel algorithm in the CHOPCHOP pipeline (see 
Labropoulos et al, submitted). In order to sample the large structure 
of the foregrounds at scales between 5 and 10 degrees we need in- 
terferometer spacing between 6.5 and 13 meters. Thus the PSF acts 
as a high-pass spatial filter. 



Figure 8 shows 'dirty' maps of the simulated Galactic syn- 
chrotron emission (MODEL B) observed with the core stations of 
the LOFAR telescope at 138 MHz. The total and polarized inten- 
sity maps are shown in Fig. 8.a & b, while the polarization an- 
gle is presented in Fig. 8.c. Note that the large scale structures of 
the emission are missing as the smallest baseline length is approx- 
imately 50 m. 



7 CONSEQUENCES FOR 21-CM REIONIZATION 
DETECTION EXPERIMENTS 

One of the major challenges of the EoR experiments is the extrac- 
tion of the EoR signal from the astrophysical foregrounds. The ex- 
traction is usually formed in total intensity along the frequency di- 
rection due to the following characteristics: 

(i) the cosmological 21 cm signal is essentially unpolarized and 
fluctuates along the frequency direction (see Fig. 7) 

(ii) the foregrounds are smooth along the frequency direction 
in total intensity and should only show fluctuations in polarized 
intensity (see Fig. 6, an example of the Galactic emission that is a 
dominant foreground component). 

Thus, the fluctuating EoR signal can be extracted from the fore- 
grounds by fitting the smooth component of the foregrounds out 
(e.g. see fig. 12 in Jelic et al. 2008). 

All current EoR radio interferometric arrays have an instru- 
mentally polarized response, which needs to be calibrated. If the 
calibration is imperfect, some part of the polarized signal is trans- 
ferred into a total intensity and vice versa (hereafter 'leakages'). As 
a result, the extraction of the EoR signal is more demanding. 

Moreover, the polarized signal could have similar frequency 
fluctuations as the cosmological signal and as such could possibly 
severely contaminate it. Thus, to reliably detect the cosmological 
signal it is essential to minimize the 'leakages' and to observe in the 
regions with the weak polarized foreground emission. We illustrate 
this through an example for the LOFAR telescope, but the problem 
is common to all current and planned EoR radio arrays. 

7.1 Levels of 'leakage' 

The 'leakages' of the total and polarized signal are produced by 
two effects: the geometry of the LOFAR array and the cross-talk 
between the two dipoles in one LOFAR antenna. The cross-talk, a 
leakage in the electronics that can cause the power from one dipole 
to be detected with other, is small compared to the geometric effects 
(Labropoulos et al, submitted) and we will ignore it for the purpose 
of this paper. 

The geometry of the LOFAR telescope is such that the array 
antennae are fixed to the ground. Therefore, the sources are tracked 
only by beam-forming and not by steering the antennae mechan- 
ically towards the desired direction. This implies that, depending 
on the position of the source on the sky, a non-orthogonal (except 
at the zenith) projection of the two orthogonal dipoles is visible by 
the source. This projection further changes as the source is tracked 
over time. Thus, the observed Stokes brightness of the source. Sobs, 
is given by (Carozzi & Woan 2009): 



Sobs = MS, 



(20) 



where M is a Mueller matrix that quantifies the distortions of a 
true source brightness S = (/, Q, [/, V) based on above geometry- 
projection effect. The Mueller matrix is defined as: 
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Figure 7. A random line of sight through a simulated 21 cm data cube for the 'Stars' patchy reionization history model (solid line). Dotted line shows the 
'leakage' of the polarized Galactic emission to the total intensity and dashed line is a sum of the two. We assume 0. 15% residual 'leakage' and we use model 
D as an example of the Galactic emission. The angular and frequency resolution of the data match that of the LOFAR telescope. 
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with the assumption of a coplanar array. Note that (/, m, n) are di- 
rection cosines and that the level of geometry-projection 'leakage' 
therefore varies across the map. 

The calculated leakages, for the LOFAR telescope observing 
at 138 MHz a 5° x 5° patch of the sky around the zenith, vary 
between 0.1% - 0.7% across the image. For the sky model we use 
the total and polarized intensity maps of Galactic emission (model 
C). Further, an instant imaging is assumed, i.e. the sky is not tracked 
over time. Note that the leakages are tiny around the center of the 
image, but they increase towards the edges. 

The same calculation we repeat for an patch of a sky at 45° 
altitude. The leakages are now much larger, e.g. for the center of 
the image the leakage to the total intensity is around 2%, but can 
reach 20% towards the edges of the image. 

Once the tracking of the sources is taken into account, the 
leakages become even more significant towards the horizon and on 
the axes parallel to the dipoles, i.e. varying between 0.1%-100% 
across the sky (for details we refer to Labropoulos et al, submit- 
ted). 

However to be as realistic as possible, we need to take into 
account that for calibration of the instrument one can use a model of 
the beam. Using the beam model, one can correct for the polarized 
instrumental response and decrease the level of 'leakages'. Thus, in 
further discussion we refer only to the residual 'leakages'. Note that 
for the purpose of this paper we assume that the residual 'leakages' 
are l%-5% over the whole field of view. 



7.2 'Leakage' of the polarized foreground 

Here we would like to point out that the residual 'leakages' caused 
by the geometry-projection effect are significant in terms of EoR 
signal detection. If these 'leakages' are not taken properly into ac- 
count during the calibration of the instrument, the polarized Galac- 
tic emission could creep into total intensity signal and severely con- 
taminate the EoR signal. This is illustrated in the Fig. 7. 

Figure 7 shows a random line of sight through a simulated 
21 cm data cube for the 'Stars' patchy reionization history model 
(solid line). As noted before, the EoR signal shows fluctuations 
in total intensity. Dotted line shows the 'leakage' of the polarized 



Galactic emission in the total intensity and dashed line is a sum of 
the two. 

Since we have chosen the Galactic emission model with the 
differential Faraday rotation (model C), the 'leakage' of the Galac- 
tic emission shows structures along the frequency direction similar 
to the EoR signal. Note that the Galactic polarized emission is as- 
sumed to be ^ IK (in this case our model C), while the level of 
the residual 'leakage' is assumed to be 1.5%. The 'leaked' Galac- 
tic polarized emission is then of the same order as the EoR signal 
(- 15 mK). 

In other words, the right combination of the Galactic polar- 
ized emission and an inaccurate calibration can result in structures 
in total intensity that have similar characteristics of the EoR signal. 
Without knowing the exact characteristics of the Galactic polarized 
emission in that region it would be impossible to extract the EoR 
signal. Therefore, the observational windows for the EoR experi- 
ments need be in regions of the Galaxy that have very weak or no 
polarization. Assuming the same residual 'leakage' of 1.5%, the 
Galactic polarized emission should not be stronger than ^ 0.1 K. 
The 'leaked' foreground contamination would then be at least an 
order of magnitude lower than the EoR signal, which is acceptable 
for the EoR detection. 

However, note that this upper limit on the polarized contami- 
nation highly depends on morphology of the Galactic emission. For 
example, if the Galactic polarized emission does not show struc- 
tures along the frequency direction (like in our model A), then any 
'leaked' contamination can be fitted out in the same way as the 
foregrounds in total intensity. Moreover, if the Galactic polarized 
emission shows significant structures on scales smaller then reso- 
lution of the array, beam depolarization will lower the level of the 
observed polarized foregrounds and therefore lower the 'leaked' 
contamination as well. 

Other possible ways of eliminating a 'leaked' polarized fore- 
ground from the data are: (i) identifying the polarized emission in 
rotation measure space using a Faraday rotation measure synthesis 
(Brentjens & de Bruyn 2005); (ii) using polarization surveys ob- 
tained by a different radio telescope and with angular resolution 
that is preferably higher than the one of the EoR experiments; and 
(iii) observing the EoR in a multiple regions of the sky. 

Note that if someone uses a different polarization data as 
an EoR foreground template, an instrumental beam depolarization 
should be take into account. At the same time the advantage is that 
radio telescopes, which are used for the polarization surveys, suffer 
from a different systematics (including 'leakages') and therefore 
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Figure 8. 'Dirty' maps of the simulated Galactic synchrotron emission (model B) observed with the core stations of the LOFAR telescope. The total and 
polarized intensity maps are shown in Fig. 8.a and Fig. 8.b, while the polarization angle is presented in Fig. 8.c. The images are simulated at 138 MHz 



comparison of the data observed with a different radio telescopes 
helps in analysis. 

The observations of the EoR signal in multiple windows have 
advantage of allowing us to cross-correlate the data and increase 
the significance of the EoR detection. Since the Galactic emission 
varies across the sky, each observational window will have its own 
characteristic foreground emission and 'leakages'. Residuals after 
foreground removal should not correlate between different obser- 
vational windows, while the EoR properties should be similar. 



8 SUMMARY AND CONCLUSIONS 

This paper presents Galactic foreground simulations used as tem- 
plates for the LOFAR-EoR testing pipeline. The simulations pro- 
vide maps of the Galactic free-free emission and the Galactic syn- 
chrotron emission both in total and polarized intensity. The maps 
are 10° x 10° in size, with ^ 1 arcmin resolution and cover the fre- 
quency range between 115 and 180 MHz pertaining to the LOFAR- 
EoR experiment. The code however is flexible as can provide sim- 
ulation over any scale with any spatial and frequency resolution. 

The Galactic emission is calculated from a 3D distribution of 
cosmic ray and thermal electrons, and the Galactic magnetic field. 
The model assumes two magnetic field components: regular and 
random. The latter magnetic field and the thermal electron den- 
sity are simulated as Gaussian random fields with power law power 
spectra. In addition, the spatial variations of the energy spectral in- 
dex p of the cosmic ray electrons are introduced to mimic the ob- 
served fluctuations of the brightness temperature spectral index (3. 
Note that all parameters of the simulation can be tuned to any de- 
sired value and this allows to explore the whole parameter space. 

The total and polarized Galactic maps are obtained for four 
different models of Galactic emission (see Fig. 3 & 4). The first 
assumes that synchrotron and free-free emitters are spatially sep- 
arated, such that thermal plasma acts as a "Faraday screen". The 
amplitude of the polarized emission is unchanged, while the polar- 
ization angles Faraday rotate. Other three simulation have regions 
where both types of emitters are mixed in different ways. The syn- 



chrotron emission is differentially Faraday rotated and depolariza- 
tion occurs (see Table 1). 

The main result of our simulations is that we are able to pro- 
duce realistic Galactic polarized emission that is comparable to ob- 
servations, i.e., presence of the structures at different scales, spa- 
tial and frequency variations of the brightness temperature and its 
spectral index, complex Faraday structures, and depolarization. The 
importance of this result comes from the fact that the planned EoR 
radio arrays have a polarized response and the extraction of the 
EoR signal from the foregrounds is usually performed along the 
frequency direction. The Galactic foreground is a smooth function 
of frequency in a total intensity and it can show fluctuations in 
polarized intensity. The EoR signal is expected to be unpolarized 
and to show fluctuations along the frequency direction. Therefore, 
an imperfect calibration of the instrumental polarized response can 
transfer a fraction of the polarized signal into a total intensity. As 
a result, the leaked polarized emission can mimic the cosmological 
signal and make its extraction almost impossible (see Fig. 7). 

Based on our simulations, we conclude that the EoR obser- 
vational windows need to be in regions with a very low polarized 
foreground emission, in order to minimize 'leaked' foregrounds. 
Faraday rotation measure synthesis, polarization surveys obtained 
by different radio telescopes, and a multiple EoR observations will 
help in mitigating the polarization leakages. However, further simu- 
lations and observations are necessary to pin point the best strategy 
for the EoR detection. 
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